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We study the nuclear spin-lattice relaxation rate l/Ti in the two-leg antiferromagnetic spin-1/2 
Heisenberg ladder. More specifically, we consider the contribution to l/Ti from the processes with 
momentum transfer (tt, tt). In the limit of weak coupling between the two chains, this contribution is 
of activation type with gap 2A at low temperatures (A is the spin gap), but crosses over to a slowly- 
decaying temperature dependence at the crossover temperature T « A. This crossover possibly 
explains the recent high-temperature NMR results on ladder-containing cuprates by T. Imai et al. 
QQ \ Recent NMR experiments on CU2O3 ladders in A14CU24O41 compounds (Ai4=La6Ca8, Sri4, SrnCaa) reveal un- 
0^ ' expected behavior of the nuclear spin-lattice relaxation rate 1 /Ti at temperatures of order of the spin gap 0] . While 
On ] both and exhibit activation-type behavior at low temperature, a well-pronounced crossover is observed 

in ^'^Cu nuclear spin relaxation rate. There exist numerous theoretical studies of 1/Ti in antiferromagnetic two-leg 
ladders , but none of them predicts a similar crossover. In this paper we restrict our attention to the undoped 
spin ladder and argue that the observed crossover is due to the processes with the momentum transfer q — {tt,tt). 
[ We describe the spin ladder by the Heisenberg Hamiltonian 
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C ' Here n is an integer labeling the rungs of the ladder, the superscripts refer to the two chains. J and J± are the 
. coupling constants. From analyzing experimental data on magnetic succeptibility Q and on ^^O NMR Knight shift 

it has been suggested that in CU2O3 ladders Jj_/J ~ 0.5. 
. Both the strong-coupling {J±/J S> 1) and the weak-coupling {J±/J ^ 1) approaches confirm that the low-lying 
00 ' magnetic excitations in the ladder are spin-1 magnons with the minimal gap at fc = (tTjTt); in the physically relevant 
On range J±/J ~ 0.5, the value of the gap is known to be A sa 0.5Jj_ We beheve that the weak-coupling 

limit captures the characteristic features of the nuclear spin relaxation. For the sake of simplicity we assume that 
T, J_L <C J (with no assumptions on the relative magnitudes of J_l and the temperature T) . Then all excited magnons 
have momenta close to the wavevector k = (tt^tt). The relaxation rate l/Ti is given by 

where F{q) are the appropriate coupling constants, loq is the nuclear resonance frequency (wq ^ ^^T), S{q,uj) is the 
dynamical structure factor, and x(g, uj) is the dynamical magnetic succeptibility. We employ the system of units with 
h = kg = = 1. Since the excited magnons have momenta close to (tt, tt), there are two major contributions to 
■ l/Ti: that with q « (0,0) (via even-magnon-number processes) and that with q « (7r,7r) (via odd-magnon-number 
processes). Accordingly, define 

(defined this way, the quantities (l/Ti)q^o ^-nd (1/Ji)q^7r do not have the dimension of inverse time). Then the total 
relaxation rate 1/Ti is a linear combination of (l/Ti)q^o ^^^id (1/T'i)g=7r. 
The coupling constants F{q) for the oxygen sites vanish a.t q ~ {tt, tt) [Q, so 

The experiments on SrCu203 suggest that in CU2O3 ladders the coupling constant F{q) for the copper site is dominated 
by the on-site hyperfine interaction and therefore only weakly depends on g ||^. Thus we expect that is a linear 

combination of (l/ri)g=o and (l/ri)g=7r with the coefficients of the same order of magnitude. Next, the available 
numerical studies M indicate that at temperatures of order of the spin gap the contributions {l/Ti)q=o and {l/Ti)q=^ 
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have comparable magnitudes. Therefore, we suggest that in the experimentally relevant temperature range is 

not dominated by the q = contribution, in contrast with other studies |^,^. 

A simple argument shows that the contribution (1/Ti)g=7r has gap 2A as opposed to gap A for (l/Ti)g=o (this 
argument was originally proposed for spin-1 chains JTo| , but it also remains valid for spin-1/2 ladders). In (1/Ti)q=0; 
the nuclear spin is relaxed by the quasi-elastic scattering of a thermally excited magnon, which requires a gap of A. 
On the other hand, (l/ri)^^^ assumes an elastic scattering with momentum transfer near (tt, tt), which means that 
the total number of incoming and outcoming magnons must be odd (each magnon carries momentum close to (tt, tt)). 
Processes of creating and annihilating a single magnon cannot occur at zero energy because of the gap. Therefore, 
the leading contribution comes from three-magnon processes: two thermally excited magnons are converted into a 
single magnon carrying the total energy of the incoming magnons (up to the NMR frequency loq), or vice versa (see 
Fig. 1 below). Such processes require energy at least 2 A. 

On the basis of the above argument, the contribution (1/Ti)q=7r was usually neglected. However, several recent 
works conclude that the effective gap in (l/Ti)g=o is somewhat larger than A (determined from succeptibility), 
either because of magnon interaction |^ or due to the singlet excitation mode |^. A larger gap in 1/Ti than in the 
succeptibility is indeed reported in CU2O3 ladders |11 |]. In view of these results, our suggestion of importance of 



q = TT contribution to ^'^Cu nuclear spin relaxation rate appears more plausible. 

For the rest of the paper we focus on computing (1/Ti)g=7r in the weak-coupling limit J± <C J. Although we 
are unable to find a closed analytic form for {l/Ti)q=T^ in the whole range of temperatures, we find the high- and 
low-temperature asymptotics and estimate the crossover temperature. 

According to the results of in the weak-coupling limit, the spin-1/2 two-leg ladder is equivalent to four massive 
Majorana fermions, combined into a triplet of mass A and a singlet of mass 3A: 

Hf^HA [6 ]+Ha[C2]+ Ha [6] + H3A [p] + , (5) 

where each H„i[^] is a free massive Hamiltonian of a Majorana fermion ^; Hint is the four-fcrmion interaction arising 
from the marginal term in the interchain coupling It has been argued in |^ that the interaction would lead to 
nonperturbative effects in (l/Ti)q=o- However, the staggered magnetization is nonlocal in terms of fermions, and we 
expect that the interaction will play a smaller role in (l/Ti),^^. Following we neglect Hint- 

In the continuum limit, the local magnetization 5*^"-' may be split into the uniform and staggered components: 

= J^") + (-!)"«("). (6) 

While J has a quadratic expression in terms of the Majorana fermions, n is non-local in the fermion operators. If the 
Majorana fermions are mapped onto non-critical l-|-l-dimensional Ising models (we need four Ising models — one for 
each Majorana fermion), the operator n may be expressed in terms of order (a) and disorder (/i) parameters of the 
Ising models j8|. In particular, 

n~ = n'P — n{/^^ cx aia2H3cr*, (7) 

where a* is the order parameter for the Ising model with gap A* = 3A; cti, (72, and /is — order and disorder 
parameters of the three identical Ising models with gap A [|j . Conventionally, we have chosen the Ising models to be 
in ordered state, so that (a) ^ at zero temperature. 

The expression (j^) would enable us to compute {l/Ti)q^Tr using (§, should we know the non-critical Ising correlation 
functions at finite temperature. In a recent work p^ , the latter has been expressed as a series in number of fermionic 
excitations, and we shall use their result to match the high- and low-temperature asymptotics of the Ising correlation 
functions. 

The prportionality coefficient in (^ is non-universal. To fix the relative normalization of the high- and low- 
temperature asymptotics, define 

= CTiCT2/i3Cr* (8) 

with the order and disorder operators normalized by their short-distance asymptotics: 

(a(x)a(O)) ^ {fi{xMO)) ^ x ^ 0. (9) 

In what follows we normalize (l/ri)g=7r accordingly, in other words, in (^) we set S{x,t) = {nj {x,t)nj (0,0)) . 
In a natural way, the following three limits may be distinguished: 
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(i) T ^ A. In this low-temperature limit, {1/Ti)q=-^ has activation-type behavior with gap 2A. The prefactor may 
be computed using the results of [ p^ . 

(ii) T 3> A*. In this limit the gaps presumably play no role, and the system is equivalent to two uncoupled 
chains ("quantum critical" phase). The result for (l/ri)g=^ in this limit may be borrowed from |13|. Although this 
high temperatures are beyond the experimental range, considering this limit is useful for determining the degree of 
applicability of the approximations being made. 

(iii) A <C T <C A*. In the ladder system this intermediate limit is never attained, since A* = 3 A. However we may 
formally take this limit assuming A* 3> A. The validity of this approximation will be discussed further. 

Below we describe in more detail the computations and the results in these three limits. 

(i) T <C A. According to the (non-ordered) Ising correlation functions may be expressed as the series in the 
number of fermionic excitations: 

i/ run ^TT 

n oven Si i—1 

where = ±1 (z — 1, . . . , n) are the particle- hole indices, 

Ml3) = [^ + oxp{-eE{l3)/T)]-' (11) 
is the Fermi distribution function, Pi are the rapidities of the excitations with energies £{(3) = A cosh f3. 

I LcliliJ 

i<j 



F(/3i, . . . , /?„)e„...,e„ = n tanh (12) 



are the corresponding formfactors p^ , p^ . 

The same expression gives (//(t, 0)/i(0, 0)), with the only difference that the sum is taken over odd numbers of 
excitations n. 

The zero-temperature magnetization ctq may be taken from the exact result on the Ising model |jl^,|l6| : 

a2 = Ai/42i/6g-i/4^3 ^ 1.8437A1/4, (13) 

where A = exp[l/12 - C'(-l)] = 1.282427 ... is the Glaisher constant. In the Ising model with gap 3A, cr^^ = 3i/''cr^. 
Using (noh, we express 



1 

^1 J q=7r 



+ 00 

dt{rX{x,t)rX (0,0)) (14) 



as a sum of the odd-magnon-number processes with zero energy transfer. The one-magnon process has a gap and 
does not contribute to (jlj). In Fig. 1 we show the three types of three-magnon processes contributing to (l/Ti)^^^. 

"^'^^-^ ^^^^^ 






1(2) ^ 1(2) 

(a) (b) (c) 

FIG. 1. Three-magnon processes contributing to (l/ri)q=^. 

In the diagrams, the arrows correspond to the particle-hole index Ei in ([l0|), the labels 1,2,3 — to the three different 
Majorana fermions (equivalently, to the spin of magnons). The contribution from the diagrams shown in Fig. 1 must 
be multiplied by two to account for the diagrams with arrows reversed. We neglect processes including the singlet 
channel (with gap 3A), since they have a larger gap 4A. The processes (a) and (b) in Fig. 1 can be shown to have 
temperature dependence exp(— 2A/r), while the process (c) dominates at low temperatures with the temperature 
dependence Texp(— 2A/r). Including only the latter contribution, we find: 

tX__^ - in ^ ^ ^ 8 cosh ^ cosh ^ eosh ^ " ^''^ 
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The low-T asymptotics of this expression is 



Putting in the numbers, 



^1 / q=-K 



i^Z(.W\(T\^-^_ ,16) 



(with the normaUzation (l/Ti)^^^ is dimensionless) . 

(ii) In the high-temperature hmit T 3> A*, the quantity (1/Ti)^=^ is determined by the short-distance asymptotics 
of the correlation function 

{h-{x)fi-m--. (18) 

X 

In this case Sachdev's result ||l^ gives 

(1/Ti),=, = ^ « 3.1416. (19) 

(iii) In this limit, we replace the operator a* by its zero-temperature expectation value, while pretending that the 
three remaining operators in (H) are massless. Then 

{h-{x)n-m^^^. (20) 
To handle this case, we redo the calculation of [ p^ for arbitrary exponent 77 in 

(n;(a;)n;(0))~^. (21) 



The magnetic succeptibility is then given by [T^ 



k) = — - ^ U^t)^ U UnTl (22) 



which leads to 



^ = lim / |^filmx(c.,fc) = DL^i2nTy 



(23) 



In our case (77 = 3/4, D — 



Numerically, 



2\ 



TJ,^. " r(|) 



^LA»l(2nTr'/\ (24) 



i\ 



= 7.028 - ) . (25) 

1 / q—7: \ / 



The three temperature dependences (|l7|), (p^, and ( p5| ) are plotted in Fig. 2. 

The mismatch of the asymptotics (ii) and (iii) at T w 3A is due to the roughness of the approximation (cr*(i)cr*(0)) — 
ctq^ made in (pO|). At T w A*, this correlation function decays at time scale of order T^^ (the same as the triplet- 
channel correlation functions), which substantially decreases (1/Ti)^=^. Although for A* — 3A the limit (iii) is never 
realized, it clearly indicates that at high temperature (l/ri)g=Tr slowly decreases with temperature, approaching the 
constant value of the asymptotics (ii). 
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FIG. 2. The three hmiting asymptotics of the temperature dependence of (l/ri)5=,r. Dashed Une is the quaUtative interpo- 
lation of the actual temerature dependence. 



Another qualitative consequence of our discussion is the crossover from increasing activation behavior of asymptotics 
(i) to slowly decreasing high-temperature asymptotics (iii)-(ii). The crossover occurs at T « A, in spite of the gap 
2 A of the low-temperature asymptotics. This agrees with the experimental results of T. Imai et al., who observed 
a sharp corssover in at T sa 425K in the undoped compound La6Ca8Cu2404i 0|. Moreover, the form of 

the temperature dependence of in the hole-doped compounds Sri4Cu2404i and SriiCa3Cu2404i appear 

qualitatively very close to the proposed form of (l/ri)g=^ (including a slowly-decreasing high-T behavior). We may 
speculate that this can possibly be explained by a different hyperfine coupling of the ^■^Cu spin in these compounds 
making (l/ri)g=7r contribution dominate in On the other hand, hole doping may modify the temperature 

dependence of (l/Ti)^^^, which deserves further theoretical study. 

It is worth mentioning that our interpolation of (l/ri)g=^ (dashed line in Fig. 2) is higher than the asympthotics 
(i) at low temperature. This comes from the fact that all terms in the low-temperature expansion (^0|) are positive, 
and including only the leading term underestimates (l/ri)g=,r at low temperatures. This observation ensures that 
the crossover from (i) to (iii)-(ii) is sufficiently sharp. 

To experimentally separate (l/ri)q=7r and (l/Ti)q=o contributions in the total spin-lattice relaxation rate ®'^(1/Ti), 
one can use the fact that q = n contribution does not have a singular dependence on the external field. While (l/Ti)q=o 
singularly diverges at low magnetic fields (as logiJ in the free-magnon model H or as in the spin-diffusion 

model (I/Ti)q=7r only weakly depends on the magnitude of the magnetic field. 

We hope that further experimental and numerical works will provide a better understanding of the crossover in 
(I/Ti)q^7r discussed in this paper. 

The authors would like to thank T. Imai for sharing the experimental data prior to publication, A. M. Tsvelik for 
drawing our attention to ref. |p^ , and D. H. Kim for many helpful discussions. 
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